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基于 C5G7-MOX 的 中 子 输 运 计算 方法 比较 及 
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摘要 ”特征 线 方法 (Method of Characteristics, MOC) 因 其 具备 强大 的 几何 处 理 能 力 ， 且 在 计算 过 程 中 亦 能 兼顾 
计算 成 本 和 计算 精度 ， 被 广泛 应 用 于 高 保 真 数值 模拟 计算 中 。 常 见 的 输 运 计算 方法 除 MOC 外 ,还 包括 碰撞 概率 
ik (Collision Probability Method, CP) 和 界面 流 法 〈Interface Current Method, IC) 等 。 本 文 从 方法 理论 以 及 数 
值 计算 上 将 MOC, CP 和 IC 进行 比较 分 析 ， 对 于 评估 其 在 pin-by-pin 计算 中 的 能 力 尤 为 重要 。 同 时 在 MOC it 
算 中 , 不 同 的 参数 选择 会 对 计算 成 本 和 计算 精度 产生 影响 , 因此 有 必要 进行 敏感 性 分 析 以 寻求 最 佳 参数 。 本 文江 
先 将 三 种 计算 方法 从 原理 上 进行 比较 分 析 ， 再 基于 2D C5G7-MOX 基准 题 完 成 了 数值 计算 及 MOC 参数 敏感 性 
初步 分 析 。 计 算 结 果 表 明 : MOC 在 计算 精度 、 计 算 效 率 和 内 存 开 销 上 均 优 于 CP 和 IC MOC 的 计算 耗 时 和 内 
存 开销 分 别 为 23.9min 和 37.5MB， 与 参考 解 的 相对 误差 仅 为 6.04x104。 而 CP 和 IC 的 计算 耗 时 分 别 为 MOC 的 
56.7 倍 和 15.6 倍 ， 内 存 开销 分 别 为 MOC 的 407.7 倍 和 32.8 倍 。 进 一 步 通过 参数 敏感 性 分 析 发 现 ， 网 格 划 分 对 
计算 内 存 开 销 以 及 计算 时 间 的 影响 最 大 , 极 角 的 选择 对 计算 精度 的 影响 最 大 , 并且 给 出 一 组 建议 参数 : 网 格 划分 
6x6， 极 角 为 GAUS 且 数 目 为 2， 方位 角 个 数 为 30。 该 组 参数 的 计算 耗 时 为 4$5.4min， 内 存 开销 为 264.7MB， 相 
对 误差 为 $.9x105， 归 一 化 后 的 栅 元 均 方 根 误 差 为 0.00255 。 
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Comparison of neutron transport calculation methods based on C5G7-MOX and preliminary analysis of sensitivity 
of MOC parameters 
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Abstract [Background]: Method of Characteristics (MOC) is widely used in high-fidelity numerical simulation 
because of its powerful geometric processing ability and its ability to take into account calculation cost and accuracy. 


In addition to MOC, the common transport calculation methods also include Collision Probability Method (CP) and 
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Interface Current Method (IC). [Purpose]: In this paper, the comparative analysis of MOC, CP and IC from method 
theory and numerical calculation is particularly important for evaluating their ability in pin-by-pin calculation. At the 
same time, in MOC calculation, different parameter selection will affect the calculation cost and calculation accuracy, 
so it is necessary to carry out sensitivity analysis to find the best parameter. [Methods]: In this paper, the three 
calculation methods are compared from the principle, and then the numerical calculation and the preliminary analysis 
of the sensitivity of MOC parameters are completed based on the 2D C5G7-MOX reference problem. [Results]: The 
calculation time and memory cost of MOC are 23.9min and 37.5MB respectively, and the relative error between MOC 
and reference solution is only 6.04x10^*. The computing time of CP and IC is 56.7 times and 15.6 times that of MOC, 
and the memory cost is 407.7 times and 32.8 times that of MOC, respectively. Based on the sensitivity analysis of MOC 
parameters, a set of suggested parameters are given: grid division is 6x6, pole Angle is GAUS and number is 2, and 
azimuth Angle is 30. The calculation time of this set of parameters is 45.4min, the memory cost is 264.7MB, the relative 
error is 5.9x10?, and the normalized RMS error is 0.00255. [Conclusions]: The results show that MOC is superior to 
CP and IC in accuracy, efficiency and memory cost. Through the parameter sensitivity analysis, it is found that the grid 
division has the greatest influence on the calculation memory cost and calculation time, and the choice of polar Angle 


has the greatest influence on the calculation accuracy. 


Key words Method of Characteristics, 2D C5G7-MOX, comparison of methods, sensitivity analysis, DRAGON 
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随 着 新 型 反应 堆 的 发 展 ， 更 多 非 均 匀 性 强 、 能 谱 异 构 更 加 复杂 的 堆 型 出 现 ， 为 堆 蕊 物理 计算 与 分 析 带 
来 了 挑战 。 在 高 保 真 反应 堆 物 理 计算 中 ， 常 用 的 中 子 输 运 方程 求解 方法 包括 : 特征 线 方法 (Method of 
Characteristics: MOC) ~, fil fit EAS 7 C Collision Probability Method, CP) 和 界面 流 法 (Interface Current Method, 
IC) 5$. CP 和 IC 方法 都 属于 积分 输 运 方法 ， 该 方法 最 早 由 A.Amouyal. P.Benoist 和 J. Horowitz 应 用 于 计 
算 栅 元 的 热 中 子 利用 系数 2， 并 且 得 到 了 很 高 的 计算 精度 。 因 此 积分 输 运 方法 曾 在 相当 一 段 时 间 内 被 作为 
反应 堆 物 理 计 算 中 的 一 种 重要 方法 ， 特 别 是 在 轻 水 堆 的 栅 元 能 谱 计算 以 及 燃料 组 件 能 谱 计 算 中 被 广泛 应 用 
Bl。 当 下 积分 输 运 方法 作为 中 子 输 运 求解 模块 已 在 多 个 程序 中 得 以 实现 ， 其 中 包括 APOLLO, WIMSD- 
5B6、MOSRAG 等 。 

特征 线 方法 最 早 由 英国 学 者 AskewD] 应 用 于 求解 复杂 的 输 运 计算 ， 随 着 近 些 年 来 计算 机 性 能 的 不 断 攀 
升 ， 特 征 线 方法 凭借 其 精确 的 数值 计算 和 灵活 的 几何 特性 ， 在 高 保 真 反 应 堆 物 理 分 析 中 得 到 了 广泛 的 应 用 
I。 目前 国内 外 许多 程序 都 采用 特征 线 方法 作为 输 运 求解 器 ,包括 DRAGON[%10、CASMOILI2、WIMSD3] 
等 。 特 征 线 方法 (Method of Characteristics, MOC) 是 沿 中 子 轨迹 以 特征 线形 式 求解 中 子 输 运 方程 的 一 种 方 
法 00， 该 方法 巧妙 地 运用 了 数学 技巧 ， 将 中 子 输 运 方程 用 微分 形式 表示 ， 沿 着 中 子 运 行 的 轨迹 进行 积分 求 
解 05。 特 征 线 方法 具备 强大 的 几何 处 理 能 力 ， 能 将 角度 变量 和 空间 变量 同时 进行 处 理 ， 理 论 上 能 处 理 任 意 
几何 ， 并 且 由 于 特征 线 中 线 与 线 之 间 的 独立 性 ， 使 得 特征 线 方法 具备 天 然 的 并 行 能 力 ， 在 计算 过 程 中 亦 能 
顾 计 算 效率 和 计算 精度 。MOC、CP、IC 三 种 方法 各 具 特 点 ， 但 在 进行 高 保 真 计算 时 ， 都 会 面临 计算 成 本 
的 挑战 ， 因 此 有 必要 对 以 上 三 种 输 运 计算 方法 进行 分 析 比 较 ， 评 估 其 在 pin-by-pin 计算 中 的 能 

本 文 首先 比较 了 特征 线 方法 (Method of Characteristics: MOC)、 碰 撞 概率 法 (Collision Probability Method, 
CP) 和 界面 流 法 Interface Current Method, IC) 的 计算 原理 ， 并 基于 2D C5G7-MOX 基准 题 分 析 了 三 种 方 
法 的 计算 精度 、 计 算 效 率 和 内 存 开 销 。 同 时 由 于 特征 线 方法 本 身 在 计算 的 过 程 中 ， 不 同 参 数 的 组 合 会 导致 
计算 精度 、 计 算 效 率 和 内 存 开销 发 生变 化 ， 需 要 进一步 开展 MOC 参数 敏感 性 分 析 。 其 中 Mirali A 3&3 
极 角 、 方 位 角 以 及 网 格 划 分 开展 了 参数 敏感 性 分 析 ， 但 仅 研究 了 上 述 三 个 参数 对 计算 精度 的 影响 。 本 文 完 


A 


善 了 参数 敏感 性 分 析 中 的 参数 类 型 ， 对 部 分 参数 的 选取 范围 进行 了 深入 探 完 ， 分 析 比 较 了 各 种 参数 的 计算 
耗 时 、 计 算 精 度 和 内 存 开销 ， 能 够 较为 全 面 的 为 特征 线 方 法 在 计算 2D C5G7-MOX 这 类 强 非 均匀 性 组 件 问 
题 的 参数 选择 与 优化 提供 一 定 参考 。 
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ERME: AEF C5G7-MOX 的 中 子 输 运 计算 方法 比较 及 MOC 参数 敏感 性 初步 


分 析 


1 计算 方法 


1.1 特征 线 方法 


特征 线 方法 (Method of Characteristics, MOC) 的 核心 思想 是 沿 着 特征 线 方向 进行 积分 ， 将 中 子 输 运 方 


程 转化 成 特征 线 方程 ， 再 通过 对 每 一 条 平行 特征 线 进行 输 运 求解 ， 从 而 完成 对 整个 目标 区 域 的 计算 。 从 理 


论 上 讲 ， 特 征 线 方法 具有 优异 的 几何 适应 能 力 ， 只 要 特征 线 足够 密 ， 便 可 对 任意 形状 的 几何 区 域 进行 输 运 


求解 。 


使 用 特征 线 方法 进行 中 子 输 运 方程 求解 时 ， 需 将 待 求解 问题 的 几何 区 域 划分 成 众多 平 源 区 ， 并 假设 在 
所 划分 的 平 源 区 中 ， 宏 观 截 面 和 中 子 源 项 均 为 常数 。 如 图 1 所 示 ， 对 于 任意 平 源 区 六， 其 中 带 箭 头 的 平行 实 
线 分 别 为 穿 过 该 区 域 的 特征 线 ， 两 虚线 之 间 的 距离 为 第 条 特征 线 的 宽度 记 作 564mnxk， 特 征 线 沿 厦 方位 角 


为 m、 极 角 为 nx"， 其 中 截面 和 源 项 近似 分 别 如 式 CIO 和 式 〈2) 所 示 。 
YXus) 9 Xu 
Qim(s) 入 Qim 


(1) 
(2) 


nea uic du a Led 
化 的 源 项 ， 于 是 稳 态 多 群 中 子 输 运 方程 可 以 表示 为 品 


sinOn, 2 Qunnk (simnx) t NtiDisunk (simnx) Qim 


区 内 不 随 位 置 而 发 生变 


(3) 


AX O0 中 : 为 中 子 沿 着 Qmg 方 向 飞行 的 极 角 ;simmx 为 特征 线 在 x-y 平面 投影 所 得 的 长 度 ，@imnk(5) 为 


中 子 角 通 量 。 式 (3) 为 一 阶 线性 微分 方程 ， 可 以 解析 求解 得 到 下 式 : 


uk = Pin nexp(—Deisimnr) tan Lim mji mi exp(- Yeisimnk)] (4) 


X (4) 中 : p83tiix 为 出 口 处 的 中 子 角 通 d mE. XX (2) 沿 着 特征 


线段 进行 积分 ， 由 中 子 守恒 关系 可 得 : 


Pimnx Simn,k — pen $i n keXb(7XtiSimnk) TY Pun mji 7 exp( 一 Fiisimnx)|ads (5) 


P E O A A 


in out 
Qim m Pimnk— Pimn,k 


Pimnk = Yu 


XtiSimnk 


nw 


pny 
图 1 特征 线 方法 示意 图 


Fig.l Diagram of method of characteristics 


对 在 区 域 i 内 的 以 方位 角 m、 极 角 n 的 所 有 特征 线段 ， 以 体积 为 权重 进行 求 和 ， 


XXXXXX-1 


(6) 


得 其 到 区 域 i 内 的 平均 
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角 通 量 密度 ; 
Bimn = Penna 5 
R CD 中 : dua EN m ifc. m BUG SERE i KRA RERBA, VRR DOR ER 


大 小 ， 对 其 方位 角 和 极 角 进行 加 权 平 均 ， 可 得 i 区 域内 的 中 子 通 量 为 : 
pi = 5a Om DN On Pimn (8) 
式 (8) P: wm 表示 方位 角 求 积 系数 ，M 为 方位 角 的 离散 数目 ，w 为 极 角 求 积 系 数 , N 为 极 角 的 离散 数目 ， 
将 式 CD 和 式 (6) RAR (8) 中 可 得 ，i 区 域内 的 标量 中 子 通 量 密度 如 下 : 
angi | EM- Om Xni On Óima Ota PAmnk 
$i T A + Ec es nl T.I (9) 
根据 式 (9) 可 知 ， 在 计算 各 平 源 区 的 中 子 通 量 密度 时 ， 极 角 n、 方 位 角 m、 特 征 线 密度 特征 线 宽度 
54mmxk 以 及 各 平 源 区 域 划分 的 精细 程度 都 会 直接 影响 中 子 通 量 的 计算 精度 。 因 此 ， 本 文 基 于 上 述 参数 开 
展 敏感 性 分 析 ， 比 较 不 同 参数 组 合 对 计算 精度 、 计 算 时 间 和 内 存 开 销 的 影响 程度 ， 为 后 续 特 征 线 参 数 的 选 
择 和 优化 提供 参考 。 


1.2 碰撞 概率 法 

碰撞 概率 法 〈Collision Probability Method, CP) 是 求解 积分 中 子 输 运 方程 的 最 基本 方法 ， 运 用 碰撞 概 
率 法 进行 求解 时 ， 将 系统 划分 成 个 互 不 相交 的 均匀 子 区 ，V=1,2,……N， 其 中 当 ij ViVeo. BT 
Vi 划分 的 足够 小 ， 可 以 近似 地 认为 在 每 一 子 区 域内 的 截面 参数 等 于 常数 。Pi; 为 j 区 内 产生 的 一 个 各 向 同性 
中 子 ， 不 经 任何 碰撞 到 达 i 区 发 生 首 次 碰撞 的 概率 L717， 由 于 Pj 的 求解 与 j 区 的 源 中 子 或 中 子 通 量 密度 的 分 
布 有 关 ， 而 中 子 通 量 密度 正 是 所 要 求 的 未 知 量 。 因 此 需要 引入 平 源 近 似 ， 即 当 子 区 域 划分 得 较 小 时 ， 近 似 
也 认为 每 个 区 域 的 源 中 子 或 中 子 通 量 密度 的 分 布 等 于 常数 ， 同 时 为 了 简化 计算 ， 亦 要 假设 入 射 中 子 在 子 区 
表面 上 的 空间 分 布 是 均匀 的 ， 且 入 射 中 子 的 角 分 布 各 向 同性 09。 故 在 多 群 近似 的 条 件 下 ， 便 可 对 局 ) 以 及 
通 量 进行 求解 如 下 : 


x 


PLE efre -7)] i 
Pus mp de A T a Y (10) 
— 51I G E 
BeigbigV = Xj [g's1 ggpig' + Sio|PuaV; ,9 = 12G (11) 


X (10) 中 : Pijg 为 g 群 内 在 j 区 内 产生 的 各 向 同性 的 中 子 ， 不 经 任何 碰撞 到 达 i 区 发 生 首次 碰撞 的 概率 。 
Fig 为 8 群 中 子 鹤 面 近似 得 到 的 宏观 截面 ，r” 一 + 为 空间 上 位 置 的 变化 ，tg (r” 一 + ) 为 从 r' 到 7 点 之 间 的 
自由 程 ，V 为 区 域 j 的 体积 。 式 (11) 称 为 碰撞 概率 形式 的 中 子 积分 输 运 方程 ， 其 中 等 号 右 端 方 括号 内 项 表 
示 在 态 内 单位 体积 内 g 群 的 源 中 子 数 ， 因 此 ， 等 号 右 端 整体 便 表示 所 有 7 区 的 g 群 中 子 在 Vi 区 域内 首次 发 生 
碰撞 的 数目 ， 等 于 等 号 左 端 丰 内 的 碰撞 数 品 ;ygis 灰 。 且 在 整个 碰撞 概率 法 求解 的 过 程 中 ， 中 子 飞 行 首次 磁 
撞 概 率 和 矩阵 忆 的 求解 最 为 复杂 和 耗 时 品 。 


1.3 界面 流 法 

界面 流 法 (Interface Current Method, IC) 的 基本 思想 是 : 将 系统 分 为 若干 子 区 ， 对 每 个 子 区 进行 均匀 
化 处 理 后 ， 再 对 每 个 子 区 建立 中 子 输 运 方程 ， 各 子 区 之 间 通 过 界面 上 的 中 子 流连 续 的 条 件 进 行 耦合 ， 以 此 
来 求解 中 子 通 量 。 

将 系统 划分 为 NIXNJ 个 子 区 ， 假 定 中 子 源 的 角 分 布 为 各 向 同性 ， 空 间 分 布 采 用 线性 分 布 ， 即 所 谓 的 线 
性 源 近似 ， 于 是 中 子 源 项 可 表示 为 07; 


y M 2 G| +a (2-2) +a (2-2, 0x «0v «o (12) 


XX C2) P: Q(x,y,0) 为 子 区 内 点 (x,y) 人 处 沿 0 方 向 的 中 子 源 密度 ，Q 为 子 区 的 平 源 强 ，a、b 为 子 区 在 x、y 
JARAK, ap ay HME x y 方向 上 的 线性 分 布 系数 。 


一 
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如 图 2 的 矩形 几何 为 例 ， 对 于 子 区 (i, 让 设 体积 为 Vi;， 根 据 中 子平 衡 可 得 在 界面 流 方程 为 : 
ZKk=4BCD zz Oa Jia — X42 Jia) + «goij = Qu (13) 


X 030 P: 9ij 为 子 区 G, D RFID, LATE G D 的 总 截面 ，Qi T IX G, D 内 的 平均 源 
强 ( 包 括 散 射 源 、 和 裂变 源 、 外 中 子 源 )。E 为 矩形 几何 子 区 的 四 个 界面 CA, B, C, DO 如 图 2 所 示 ，9 为 各 
界面 上 中 子 飞 行 方向 划分 为 四 个 象限 ， 每 个 边界 的 四 个 象限 在 图 2 中 均 己 标 出 ， 厌 ,为 子 区 界面 上 象限 q 的 
出 射 中 子 流 ， Jea FX AH k 象限 g 的 入 射 中 子 流 。 

界面 流 法 在 求解 通 量 时 ， 考 虑 的 是 每 个 子 区 域内 源 与 相 邻 子 区 域 的 边界 的 中 子 守恒 关系 271， 每 个 子 区 
只 需 计 算 首 次 穿 透 和 泄露 两 个 概率 ， 且 计算 的 概率 数目 只 与 系统 内 所 出 现 的 子 区 类 型 成 正比 。 而 碰撞 概率 
法 中 首次 磁 撞 的 计算 量 ， 却 是 与 系统 内 子 区 数目 的 平方 成 正比 。 因 此 界面 流 法 相 较 于 碰撞 概率 法 ， 具 有 较 
高 的 计算 效率 。 


LI 


图 2 象限 的 划分 与 界面 标号 


Fig.2 Quadrant division and interface labeling 


2 计算 模型 


2.1 2D C5G7 MOX 基准 题 
Reflection Boundary 


Reflection Boundary 
Vacuum Boundary 


Vacuum Boundary 
3 C5G7-MOX 组 件 布置 示意 图 
Fig.3 C5G7-MOX assembly layout diagram 


C5G7 基准 题 是 由 OECD/NEA 于 2003 年 发 布 的 基准 题 ， 主 要 用 于 评价 现 有 输 运 计算 方法 进行 全 堆 芯 
非 均 匀 计 算 的 能 力 09。 该 基准 题 包 含 y UOM MOX 两 种 燃料 组 件 ， 每 个 组 件 的 边 长 为 21.42cm， 由 17x17 
的 燃料 栅 元 组 成 。 外 围 慢 化 剂 区 域 的 宽度 与 组 件 边 长 相同 。2D C5G7-MOX 基准 题 的 结构 如 图 3 Hz. dh 
侧 和 下 方 为 真空 边界 ， 左 侧 和 上 方 为 全 反射 边界 19。 各 类 组 件 中 的 各 个 栅 元 都 具有 相同 的 几何 参数 。 栅 元 
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几何 结构 如 图 4 所 示 , 其 中 包括 燃料 和 包 壳 打 混 的 圆 形 区 域 以 及 外 部 慢 化 剂 的 矩形 区 域 , 栅 元 边 长 为 1.26cm， 
燃料 区 域 的 半径 为 0.54cm。 

C5G7-MOX 基准 题 共 由 七 种 材料 构成 : UOo 燃料、4.3%、7.0% 和 8.7% 的 MOX 燃料 、 裂 变 室 、 导 向 管 
和 水 。 其 中 水 不 仅仅 用 于 填充 外 部 慢 化 剂 区 域 ， 还 用 于 填充 每 个 机 元 内 部 棒 外 的 区 域 ,， UOs 和 MOX 组 件 
内 的 燃料 棒 布置 如 图 5 所 示 。 


Moderator 


Fuel-Clad Mix 


1.08 
1.26 


4 C5G7-MOX 栅 元 结构 示意 图 
Fig.4 C5G7-MOX cell structure diagram 


图 5 C5G7-MOX 燃料 棒 布 置 示意 图 
Fig.S C5G7-MOX fuel rod layout diagram 


22 DRAGONS.0.8 

DRAGON 是 由 加 拿 大 蒙特 利 尔 理工 大 学 核 工程 研究 所 开发 的 , 基于 确定 论 方法 的 反应 堆 栅 格 计算 程序 
外 。 该 程序 采用 特征 线 方法 、 界 面 流 法 、 碰 撞 概 率 法 等 多 种 数值 方法 求解 中 子 输 运 方程 PZ0， 已 被 广泛 用 于 反 
应 堆 物 理 计 算 研 究 中 。DRAGON 具有 较 强 的 几何 处 理 能 力 ， 能 处 理 常 见 的 二 维和 三 维 几何 模型 。 基 于 
DRAGON 计算 2D C5G7-MOX 基准 题 的 方法 如 图 6 所 示 。 首 先 通过 MAC 模块 读 入 该 基准 题 中 各 材料 的 7 
群 截面 参数 ， 几 何 上 采用 51x51 的 栅 元 阵列 排 布 ， 其 中 顶部 和 左 侧 的 边界 条 件 为 镜面 边界 ， 右 侧 和 底部 的 
边界 条 件 为 真空 边界 〈GEO 模块 ) ， 然 后 对 待 求解 的 几何 区 域 进行 特征 线 布 置 与 追踪 ， 获 得 各 平 源 区 内 特 
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内 的 中 子 通 量 密度 及 反应 率 ， 完 成 二 维 输 运 计算 (ASM 模块 、FLU 模块 和 EDI 模块 ) 。 


Macrolib 
Geometry 


Ameliorated 
track 
Asmpij 


图 6 基于 DRAGON 的 特征 线 方法 
Fig.6 Method of characteristics based on DRAGON 


Geometry 
track 


Flux 


3 计算 结果 及 分 析 


本 文采 用 蒙特 卡 罗 程 序 OpenMC 的 计算 结果 作为 参考 解 。OpenMC 是 2011 年 由 美国 肪 省 理工 学 院 开 


征 线 的 相关 几何 信息 CNXT 模块 和 MCCGT 模块 ) ， 最 后 逐一 对 每 条 特征 线 进行 输 运 求解 ， 获 得 各 平 源 区 


TUN 


发 的 蒙特 卡 罗 ! 


子 和 光子 输 运 模拟 程序 ， 专 注 于 反应 堆 建 模 和 反应 堆 物 理 方法 研究 ， 旨 在 实现 高 保 真 


性 能 、 大 规模 仿真 计算 21223， 有 具备 临界 计算 、 固 定 源 计算 等 多 种 功能 23]。 
测试 环境 为 PC 单机 ， 有 具备 5.0GHz 十 四 核心 Intel Ef i9 处 理 器 ，16GB 内 存 。 通 过 OpenMC 共计 算 


^ [可 


模拟 800 代 ， 去 掉 前 200 代 ， 统 计 后 600 代 ， 每 代 跟 踪 600000 个 中 子 ， 以 得 到 各 燃料 栅 元 的 归 一 化 功率 及 


X (14) P: k 


D 为 本 文 计算 的 有 效 中 子 增殖 因子 ， 而 Ko 


有 效 中 子 增殖 因子 (kef 王 1.18652+0.00004)。 相 对 误差 计算 公式 如 下 : 


_ [kekse 
RE 一 e 4 e 


keff 


eff 


为 OpenMC 计算 的 有 效 中 子 增殖 因子 。 通 过 


(14) 
比较 


组 件 中 每 个 栅 元 的 功率 计算 误差 ， 从 而 进一步 对 功率 的 空间 分 布 进行 评估 ， 均 方 根 误差 计算 公式 如 下 : 
RMSE = [Eco 
N 


X (15) P: N 为 栅 元 的 总 数目 


的 归 一 化 功率 。 


3.1 不 同 中 子 输 运 方法 计算 比较 


首先 将 特 生 


E 线 方法 (MOC) 与 碰撞 概率 法 CCP) 、 界 面 流 法 AC) ， 通 过 对 2D C5G7-MOX 基准 


(15) 


，pp ;为 本 文 计算 所 得 第 ;个 栅 元 的 归 一 化 功率 ，po; 为 参考 解 中 第 ;个 栅 元 


表 1 展示 了 分 别 采用 特征 线 方法 、 碰 撞 概 率 法 和 界面 流 法 在 计算 2D C5G7-MOX 基准 题 的 耗 时 、 


开销 以 及 相对 误差 。 计 算 时 间 方 面 : 特征 线 方法 耗 时 最 短 ， 界 面 流 法 次 之 ， 碰 撞 概 率 法 最 长 。 其 中 界 
法 耗 时 是 特征 线 方法 的 15.6 倍 ， 碰 撞 概 率 法 耗 时 是 特征 线 方法 的 56.7 倍 ; 内 存 开销 方面 : 特征 线 方法 
存 开 销 远 小 于 碰撞 概率 法 和 界面 流 法 ， 其 中 碰撞 概率 法 内 存 开销 最 大 。 界 面 流 法 的 内 存 开销 是 特征 线 方法 
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题 计 
算 的 计算 精度 、 计 算 效 率 和 内 存 开 销 进 行 比较 。 再 从 三 种 方法 的 计算 原理 进行 解释 ， 能 够 更 全 面 地 对 特征 
线 方法 、 础 撞 概 率 法 、 界 面 流 法 的 计算 能 力 进 行 评 估 。 


内 存 
apo 
的 内 
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的 32.8 fi, RÉJERBEXSTERUVS 103 


E 


开销 是 特征 线 方法 的 407.7 倍 ; 


法 与 参考 解 相对 误差 高 达 3.8996x10-?。 


采用 碰撞 概率 法 计算 虽然 
若干 区 域 ， 每 个 区 域 的 中 子 通 量 密度 彼此 相互 耦合 ， 得 到 中 子 通 量 需要 求解 所 有 
这 些 都 会 使 计算 规模 变 大 ， 从 而 导致 耗 时 增多 内 存 3 


也 足够 精 丰 


但 是 会 消耗 大 量 时 


计算 精度 方面 特征 线 方法 的 计算 结果 最 为 
甫 确 ， 与 参考 解 的 相对 误差 只 有 6.04x104， 碰 撞 概 率 法 算得 的 结果 与 参考 解 相对 误差 6.14x104， 而 界面 流 


碰撞 概率 法 将 系统 分 成 


FRE PE, MAEI 


于 销 增 大 。 界 面 流 法 在 计 外 


区 域 的 首次 碰撞 概率 24。 


过 程 中 ， 由 于 中 子 流 过 各 子 


区 域 之 间 的 界面 必须 在 连续 条 件 下 耦合 ， 且 每 个 子 区 域 上 只 与 相 邻 子 区 域 耦合 ， 


只 需要 足够 密 的 特征 线 将 所 有 计 
此 可 以 通过 调整 特征 线 的 疏 密 程 度 ， 在 保证 精度 的 同时 也 能 减 4 
计算 2D C5G7-MOX 基准 题 这 类 强 非 均匀 性 的 问题 时 ， 计 算 精 度 、 训 


概率 法 和 界面 流 法 。 


-« 


x1 2D C5G7-MOX 燃料 组 件 计 算 结 果 


计算 方法 耗 时 

Computing method Elapsed time / min 

MOC 23.9 

CP 1355.6 

IC 374.1 
32MOC 参数 敏感 性 分 析 


进一步 对 影响 特征 线 方法 计生 


Tablel 2D CSG7-MOX fuel assembly calculation results 


内 存 开销 


Memory overhead / MB 


37.5 


15291.1 


1233.4 


因此 在 平 源 近 似 假设 下 计生 
强 非 均匀 性 问题 的 误差 较 大 2529, 会 使 其 误差 远大 于 特征 线 方法 和 碰撞 概率 法 ,特征 线 方 法 在 计算 过 程 中 ， 
算 区 域 进行 覆盖 ， 沿 着 特征 线 进行 积分 ， 从 而 完成 对 整个 区 域 的 计算 。 因 
计算 成 本 。 综 合 以 上 分 析 ， 特 征 线 方法 在 
| 算 效 率 和 内 存 开销 总 体 上 是 优 于 碰撞 


相对 误差 


Relative error / 10? 


60.4 
61.4 


3899.6 


结果 的 关键 参数 进行 敏感 性 分 析 。 由 式 〈9) 可 知 极 角 、 方 位 角 、 特 征 线 


密度 的 选择 都 会 对 中 子 通 量 的 计算 结果 产生 影响 ， 且 每 个 栅 元 划分 的 网 格 数目 也 会 对 计算 精度 、 计 算 效率 


和 内 存 开销 产生 影响 。 


R 2 Bin: 


参数 类 型 
Parameter type 
极 角 类 型 
Polar Angle type 
极 角 数目 
Number of polar angles 
方位 角 

Azimuth angles 

线 密度 


linear density / cm'! 


通过 改变 网 格 数目 、 方 位 角 、 极 角 、 特 征 线 密度 进行 计生 
数 对 计算 效率 以 及 精度 的 影响 从 大 到 小 为 : 网 格 数 
算 产 生 的 影响 大 小 顺序 ， 分 别 对 网 格 数 


， 共 得 到 44 组 数据 ， 初 步 分 析 得 出 ， 以 上 参 


目 、 极 角 、 方 位 角 、 特 征 线 密度 进行 计生 


表 2 参数 类 型 与 选择 区 间 


Table2 Parameter type and selection interval 


选择 区 间 


Selection interval 


GAUS, LCMD. CACB 


8. 12. 16. 30. 32. 64 


8. 10. 12. 15. 20 
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目 、 极 角 、 方 位 角 、 特 征 线 密度 。 因 此 ， 本 文 按照 对 计 


分 析 ， 参 数 类 型 及 区 间 如 


获 瀚 源 等 ， 基 于 C5G7-MOX 的 中 子 输 运 计算 方法 比较 及 MOC 参数 敏感 性 初步 分 析 


网 格 划分 


Mesh generation 


2x2、4x4、6x6、8x8 


3.2.1 网 格 划 分 与 极 角 


400000 


200000 


Memory Overhead / kB 


Memory Overhead / kB 
N N N 
S 
S 


—a— CACB 2 "E, 
*— CACB 3 nad i 


4— GAUS 2 


4— LCMD3 


36 
v— GAUS 3 
*— LCMD2 

—4— LCMD 3 


16 36 64 
The Number of Individual Cell Grids 


图 7 不 同 网 格 和 极 角 下 的 内 存 开销 
Fig.7 Memory overhead for different grids and polar angles 


图 7 展示 


了 选择 不 同 的 网 格 划分 以 及 不 同 极 角 进行 计算 所 得 的 内 存 开销 ， 其 


The Number of Individual Cell Grids 


FERA es ET BIET 


格 数目 ， 纵 坐标 表示 整个 计算 过 程 得 内 存 开销 《单位 kB)， 其 中 不 同 的 曲线 分 别 代表 六 种 极 角 。 从 图 7 可 
知 ， 随 着 单个 栅 元 的 网 格 划 分 数目 的 增多 ， 内 存 开 销 也 随 之 显著 增 大 。 而 相 比 之 网 格 划分 对 内 存 开 销 的 影 
响 ， 不 同 极 角 的 选择 对 内 存 开销 的 影响 几乎 可 以 忽略 不 计 。 在 网 格 划分 为 2x2 的 


情况 下 ， 六 种 极 角 分 别 计 
算 的 内 存 开 销 为 3.8x10 和 kB; 在 网 格 划分 为 4x4 的 情况 下 ,六 种 极 角 分 别 计算 的 内 存 开 销 为 1.27x105kB; 在 


网 格 划 分 为 6x6 的 情况 下 ， 六 种 极 角 分 别 计算 的 内 存 开销 为 2.71x105kB;， 在 网 格 划 分 为 8x8 的 情况 下 ， 六 
种 极 角 分 别 计算 的 内 存 开销 为 4.36x105kB。 根 据 不 同 网 格 划 分 和 极 角 下 的 内 存 开 销 可 得 ， 相 同 网 格 数目 下 
的 不 同 极 角 划分 对 内 存 的 开销 影响 仅 在 400kB 以 内 ， 而 不 同 网 格 划分 对 内 存 开销 影响 非常 明显 ， 网 格 划 分 


越 密 内 存 开销 


更 越 大 ，8x8 网 格 划 分 比 2x2 网 格 划分 内 存 开销 多 3.98xl05kB， 是 
10000 
8000 上 
A 
s pA 
2 6000 上 Z 
E f 
$ 1000 F =— CACB 2 
a xA  |-4— CACB3 
- = —v— GAUS2 
20 r- P + GAUS 3 
t 4— LCMD2 
»— LCMD 3 
0 L fi 1 fi 


1 16 36 64 
The Number of Individual Cell Grids 


图 8 不 同 网 格 和 极 角 下 的 总 耗 时 
Fig.8 Total time spent at different grids and polar angles 
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图 8 展示 了 不 同 网 格 和 极 角 下 的 耗 时 ， 其 中 横 坐 标 表示 单个 栅 元 网 格 数目 ， 纵 坐标 表示 整个 计算 过 程 


的 耗 时 (单位 s), 其 中 不 同 的 线 分 别 代 表 六 种 极 角 。 从 图 


8 可 知 , 网 格 划分 对 耗 时 的 影响 要 大 于 极 角 选择 ， 


且 随 着 网 格 划分 数目 越 多 ， 计 算 所 需 的 耗 时 明显 增多 。 在 网 格 划分 为 2x2 的 情况 下 ， 六 种 极 角 分 别 计算 的 


平均 耗 时 为 1.3x103s; 


~ 


网 格 划 分 增加 到 4x4 时 ,六 种 极 角 分 别 计算 的 平均 耗 时 也 随 之 增加 到 2.2x103s; 在 


6x6 的 更 多 网 格 划 分 情况 下 ， 六 种 极 角 分 别 计算 的 平均 耗 时 也 上 升 到 了 2.7x103s; 直至 网 格 划分 增加 到 8x8 
时 ， 六 种 极 角 分 别 计算 的 平均 耗 时 为 7.4x103s。 图 8 中 可 知 ， 网 格 数目 从 6x6 增加 到 8x8 时 ， 耗 时 的 变化 
尤为 突出 ， 六 种 不 同 的 极 角 的 耗 时 都 会 增加 至 少 4.6x103s。 


图 9 展示 了 不 同 极 角 在 不 同 网 格 数目 下 计算 所 得 的 误差 。 不 论 选择 何 
精度 都 有 显著 增长 。 其 中 GAUS 求 积 组 
2 比 3 更 加 精确 


DRZ, ii 
选择 极 角 数目 
5 内 发 生变 化 。 但 由 
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图 9 不 同 网 格 和 极 角 下 的 计算 误差 


Fig.9 Calculation error under different grids and polar angles 


图 6 和 


3.2.2 方位 角 


显 


耗 时 有 
目 


的 增长 ， 


化 。 这 是 由 


方位 角 


Azimuth angles 


8 
12 


16 


为 8 的 3 倍 。 在 内 存 开销 
差 可 知 ， 随 着 方位 角 数 目的 增多 ， 计 和 售 
于 DRAGON 5.0.8 程序 ! 


在 基于 网 格 划 分 6x6， 极 角 选 择 GAUS(2) 
分 析 。 不 同方 位 角 下 的 内 存 开销 


但 到 了 方 


o 在 网 格 划分 分 别 为 6x6 与 8x8 时 , 364 
图 7 可 知 ， 此 时 网 格 划分 8x8 


、 耗 时 以 及 计 香 
ff y 30 之 后 耗 时 几乎 不 再 
中 可 得 出 ,改变 方位 角 的 数目 内 存 开 和 
精度 也 会 愈 发 精确 ， 
方位 角 选 择 的 上 限 为 3 


要 比 其 他 求 积 组 精确 


顾 计 算 精 度 和 效率 。 


— 


种 极 角 类 型 和 个 数 ， 随 着 网 格 划 
至 少 5x104, TE GAUS 求 积 组 中 


FE GAUSC) 极 角 的 计算 精度 仅 在 5x10- 
的 内 存 开销 是 网 格 划 分 6x6 的 1.6 倍 ， 而 耗 时 更 是 
网 格 划 分 6x6 的 2.7 倍 。 根 据 不 同 网 格 划 分 和 极 角 选择 下 对 计算 精度 、 计 算 效 率 和 内 存 开 销 的 分 析 ， 
网 格 划分 6x6， 极 角 选 择 GAUS(2)， 能 较 好 的 


4E 


导出 


的 参数 条 件 下 ， 通 过 不 同 的 方位 角 个 数 完成 进一步 的 敏感 性 
误差 如 表 3 所 示 。 从 耗 时 可 得 出 ， 随 着 方位 角 数 目的 增多 ， 
增长 ， 且 方位 角 数 目 为 30 时 耗 时 约 为 方位 角 数 
省 的 变化 在 5 多 以 内 。 根 据 表 
晶 在 方位 角 数 目 为 30 之 后 ， 精 度 便 不 再 发 生变 


3 中 的 计算 误 


表 3 不 同方 位 角 下 的 计算 结果 


0， 故 综合 考虑 方位 角 的 最 优 参数 为 30。 


Table3 Calculation results under different azimuth angles 


耗 时 
Elapsed time / s 


890.06 


1345.66 


1729.97 


内 存 开销 
Memory overhead / kB 


268492 


271336 


281766 


XXXXXX-8 


相对 误差 


Relative error / 10? 


51.16 


21.83 


13.23 
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30 2671.79 271202 
32 2723.28 271056 
64 2537.84 271412 


3.2.3 特征 线 密 度 


5.90 


5.90 


5.90 


以 网 格 划分 6x6， 极 角 选 择 GAUSO), 方位 角 为 30 的 情况 下 ， 选 择 不 同 的 特征 线 密度 进行 计算 ， 得 到 


因此 随 着 特征 线 密度 的 增 


cn 


hr 


肖 和 计算 时 间 。 


内 存 开 销 、 耗 时 及 计算 误差 如 表 4 所 示 。 不 同 特征 线 密度 下 的 计算 精度 者 
线 密 度 为 12 Cem 的 精度 最 好 , 误差 仅 为 5.9x10-5。 同时 由 于 DRAGON 程序 采用 了 特征 线段 合并 技术 ， 

加 ， 内 存 不 会 出 现 明显 的 增加 R834。 虽然 合并 了 部 分 特征 线段 ， 但 仍 会 使 得 整体 过 
旦 规模 增 大 ， 从 而 导致 耗 时 增加 。 因 此 ， 特 征 线 密度 选择 12 Cem 能 使 计算 误差 最 小 ， 同 时 兼顾 内 存 开 


BHAR, EE 3x104 以内。 其 


表 4 不 同 特征 线 密度 下 的 计算 结果 


Table4 Calculation results are obtained under different characteristic linear densities 


特征 线 密度 


Characteristic linear densi 


8 

10 
12 
15 
20 


3.2.4 功率 分 布 


耗 时 内 存 占用 相对 误差 
ties / cm! Elapsed time/s Memory overhead / kB Relative error / 10? 
1807.41 271382 14.16 
2409.81 272808 10.87 
2662.86 271268 5.90 
2856.16 273492 22.50 
3487.05 282894 12.30 


特征 线 密度 12 Ccm!) 作 


基于 OpenMC 得 到 的 归 一 化 功率 分 布 如 图 10 所 示 。 同时 以 网 格 划 分 6x6、 极 角 GAUSO 方位 角 30、 


为 选 定 参数 ， 计 算 所 得 的 功率 归 一 化 分 布 如 图 11 所 示 。 由 图 11 可 知 ， 计 算 结 果 


与 参考 解 吻 合 较 好 ， 其 ， 
列 ， 与 参考 解 的 相对 误差 
其 值 为 1.61x105, 与 参考 角 


最 大 归 一 化 功率 为 2.506， 位 于 左上 和 角 UO» 燃料 组 件 中 的 第 4 fT 5 列 和 第 5 行 4 
为 0.04%。 位 于 左下 和 角 UO» 燃料 组 件 的 裂变 室 第 26 fT 26 列 的 归 一 化 功率 最 小 ， 


坚 的 的 相对 误差 为 0.32%。 通 过 式 (15 ) 计 算得 到 整个 基 Y 


仁 题 的 均 方 根 误 差 为 0.00255， 


由 此 可 知 ， 该 参数 的 计生 


结果 与 参考 解 吻 合 良 好 。 
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图 10 OpenMC 栅 元 归 一 化 分 布 
Fig.10 Normalized cell distribution in OpenMC 


2.5 
2.0 
15 
1.0 
0.5 
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图 11 DRAGON 栅 元 归 一 化 分 布 
Fig.11 Normalized cell distribution in DRAGON 


通过 改变 不 同 参数 进行 敏感 性 分 析 可 得 ， 在 计算 精度 上 : 影响 最 大 的 因素 是 极 角 的 选择 ， 其 中 GAUS 
求 积 组 要 比 其 他 求 积 组 精确 至 少 Sx104， 由 公式 〈8) 可知， 因为 不 同 求 积 组 会 使 得 角度 离散 过 程 中 的 极 角 
求 积 系 数 w, 发生 变化， 进而 影响 最 终 的 计算 结果 R99。 在 计算 成 本 上 : 网 格 划 分 对 计算 内 存 开销 以 及 计算 时 
间 的 影响 最 大 , 原因 在 于 增多 单个 栅 元 的 网 格 数目 , 会 使 整个 51x51 的 组 件 中 每 个 栅 元 的 网 格 划 分 都 增多 ， 
从 而 显著 增加 其 计算 规模 ， 最 终 影响 到 计算 所 需 的 总 时 间 以 及 内 存 开销 。 


4 结语 


本 文 在 分 析 比 较 MOC、CP 和 IC 计算 原理 的 基础 上 ， 进一步 基 于 2D C5G7-MOX 基准 题 进行 了 数值 计 
算 分 析 , 计算 结果 表明 : 特征 线 方法 在 计算 精度 、 计 算 效率 和 内 存 开销 方面 均 优 于 碰撞 概率 法 和 界面 流 法 。 
其 中 特征 线 方法 的 耗 时 为 23.9min， 只 占 界 面 流 法 的 /16。 特 征 线 方法 内 存 开销 约 37.3SMB， 仅 为 界面 流 法 
的 1/33。 而 碰撞 概率 法 的 耗 时 、 内 存 开销 则 远大 于 以 上 两 种 方法 ， 但 特征 线 方法 和 碰撞 概率 法 的 与 基准 值 
的 相对 误差 都 约 为 6.1x104， 远 小 于 界面 流 法 的 相对 误差 。 因 此 ，MOC 在 计算 2D C5G7-MOX 这 类 强 非 均 
匀 性 组 件 问题 时 ， 其 计算 精度 、 计 算 效 率 和 内 存 开 销 总 体 优 于 碰撞 概率 法 和 界面 流 法 。 

进一步 对 特征 线 方法 进行 初步 参数 敏感 性 分 析 ， 选 择 网 格 划 分 、 方 位 角 、 极 角 类 型 与 个 数 、 特 征 线 密 
度 进 行 计算 分 析 。 通 过 不 同 参数 组 合 下 的 计算 结果 比较 了 计算 精度 、 计 算 效 率 和 内 存 开 销 。 最 终 得 出 网 格 
划分 6x6、 极 角 GAUSQ). 方位 角 30、 特 征 线 密 度 12cm- 是 一 组 能 够 较 好 的 兼顾 计算 精度 和 计算 成 本 的 参 
数 ， 该 参数 组 合 下 的 计算 耗 时 为 45.4min， 内 存 开销 264.7MB ， 相 对 误差 为 5.9x105， 归 一 化 后 的 栅 元 均 方 
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根 误差 为 0.00255， 与 参考 解 吻 合 恨 好 。 由 于 参数 的 选择 与 计算 精度 、 计 算 效 率 和 内 存 开销 间 的 关系 复杂 ， 
各 自 权 重 难以 精确 量化 评估 ， 后 续 研 究 中 拟 引 入 深度 学 习 算法 对 参数 选择 及 优化 进行 深入 研究 。 


作者 贡献 声明 ”到 瀚 源 负责 建立 模型 ， 研 究 并 分 析 数 据 以 及 起 草 论文 ， 张 彬 航 负 责 设计 具体 研究 内 容 和 
方向 ， 提 出 合理 化 研究 方案 以 及 提供 理论 文 撑 ， 同 时 指导 程序 的 使 用 以 及 数据 的 采集 ， 对 论文 的 知识 性 内 
容 做 审阅 并 且 负 责 对 论文 的 修改 ， 袁 显 宝 负责 过 程 监督 以 及 成 果 验 收 ， 张 永 红 负 责 论文 审阅 ， 唐 海 波 负责 
技术 支持 。 


XXXXXX-11 


Jh 技 术 20XX,XX: XXXXXX 


参考 文献 


1 谢 仲 生 ， 邓 力 . 中 子 输 运 理论 数值 计算 方法 [M]. 西安 : 西北 工业 大 学 出 版 社 , 2005. 
XIE Zhongsheng, DENG Li. Numerical calculation method of neutron transport theory[M]. Xi'an: Northwestern Polytechnical 
University Press, 2005.[LinkOut] 


2 Amouyal A, Benoist P, Horowitz J J. New method for determining the thermal utilization factor in a unit cell[J]. Nuclear. Energy, 
1957,6: 79. 


3 BEES. 近代 反应 堆 物理 分 析 [M]. 北京 : 中 国 原子 能 出 版 社 , 2017: 93. 
Cao Liangzhi. Analysis of modern reactor physics [M]. Beijing: China Atomic Energy Press, 2017: 93. 


4 M. Coste, H. Tellier, P. Ribon, V. Raepsaet, C. Van der Gucht, New Improvements in the Self-Shielding Formalism of the Apollo-2 
Code, Technical Report, CEA Centre d'Etudes de Saclay, 1993. 

5 T. Kulikowska, WIMSD-5B: a Neutronic Code for Standard Lattice Physics Analysis, Distributed by NEA Data Bank. Saclay, France, 
1996. 

6 K. Okumura, T. Kugo, Y. Nakano et al, “Development of Modular Reactor Analysis Code System MOSRA," Preprints L46, 2008 Fall 
Meeting of the Atomic Energy Society of Japan, (2008) [in Japanese]. 

7 Askew JR. A characteristics formulation of the neutron transport equation in complicated geometries[J]. AEEW-M, 1972, 1108. 

8 got. 基于 特征 线 方法 的 2D/1D 及 2D/3D 耦合 中 子 输 运 计算 方法 研究 [D]. 西安 : 西安 交通 大 学 , 2017. 

LIANG Liang. Research on 2D/1D and 2D/3D Coupled neutron transport calculation methods based on eigenline method[D]. 


Xi'an: Xi'an Jiaotong University, 2017. 

9 Hebert A. DRAGONS and DONJONS, the contribution of Ecole Polytechnique de Montreal to the SALOME platform[J]. Annals of 
Nuclear Energy, 2016, 87: 12-20. 

10 Le Tellier R, Hebert A. On the integration scheme along a trajectory for the characteristics method[J]. Annals of Nuclear Energy, 
2006, 33 (14-15): 1260-1269. 

11 Hsiao MY, Wheeler JK, de la Hoz C. Application of CASMO-4/MICROBURN-B2 Methodology to Mixed Cores with Westinghouse 
Optima2 Fuel[J]. Nuclear Science and Engineering, 2011, 167 (3): 230-241. 

12 Rhodes J LD, Smith K CASMO-5/CASMO-5M - A fuel assembly burnup program methodology manual. In: Report SS, editor, 2008. 

13 Lindley BA, Hosking JG, Smith PJ, et al. Current status of the reactor physics code WIMS and recent developments[J]. Annals of 
Nuclear Energy, 2017, 102: 148-157. 

14 Askew J R. A characteristics formulation of the neutron transport equation in complicated geometries[R]. United Kingdom Atomic 
Energy Authority, 1972. 

15 XRRR, TR, WER, $. 基于 DRAGON FEF RR & e Tae VES IOS EIA]. 中 国 科技 论文 , 2015, 10023): 2728 - 
2733. DOI: 10.3969/j.155n.2095 -2783.2015.23.008. 


LIU Zijing, YU Tao, XIE Jinsen, et al. Neutron transport calculation methods comparative analysis based on DRAGON code[J]. 
China Sciencepaper, 2015, 10(23): 2728-2733. DOI: 10.3969Vj.issn.2095-2783.2015.23.008.[ 知 网 中 文 ][ 知 网 英文 ] 
16 Mira M, El Hajjaji O, Jai O, et al. Derivation of optimal process MOC parameters and analysis of the 2D C5G7 MOX benchmark 


using DRAGONS code[J]. Nuclear Engineering and Design, 2022, 388: 111613. DOI: 
10.1016/j.nucengdes.2021.111613.[LinkOut] 
17 XAR. 中 子 输 运 方程 确定 论 数值 方法 [M]. 北京 : 中 国 原子 能 出 版 社 , 2018: 65. 
Wu Hongchun. Numerical method for determining Neutron Transport Equation [M]. Beijing: China Atomic Energy Press, 2018:65. 
18 谢 仲 生 . 中 子 输 运 理论 数值 计算 方法 [M]. 西安 : 西安 交通 大 学 出 版 社 , 2022: 158. 
Xie Zhongsheng. Theoretical Numerical Calculation Method for Neutron Transport [M]. Xi 'an: Xi 'an Jiaotong University Press, 
202:158. 
19 Marleau G, Hebert A, Roy R. A user guide for DRAGON version 4, Report IGE-294[Z]. Montr'eal: Institut de génie nucléaire, école 
Polytechnique deMontréal, 2015 
毕 光 文 , 司 胜 义 .DRAGON 程序 用 于 压 水 堆 燃料 组 件 计 算 的 有 效 性 验证 与 分 析 [C]V 第 十 三 届 反 应 堆 数值 计算 与 粒子 输 运 学 
术 会 议 凤 2010 年 反应 堆 物 理会 议论 文集 . 西安 : 第 十 三 届 反 应 堆 数值 计算 与 粒子 输 运 学 术 会 议 , 2010: 1-8. 
BI Guangwen, SI Shengyi. Validation and analysis of calculations for PWR assemblies with lattice code DRAGON[C J//Proceeding 


T 


I 


XXXXXX-12 


RE: IT C5G7-MOX 的 中 子 输 运 计算 方法 比较 及 MOC 参数 敏感 性 初步 分 析 


sof13th Reactors and Particle Transport Numerical Conference and Reactor Physics Conference in 2010. Xi'an: 13th Reactors and 
Particle Transport Numerical Conference and Reactor Physics Conference, 2010: 1-8. (inChinese) 

21 AIEK 基于 OpenMC 的 空间 核反应 堆 临 界 及 安全 分 析 [D]. 江西 : 东 华 理工 大 学 , 2022. 
Hao Pengfei: Criticality and Safety Analysis of Space nuclear reactor based on OpenMC [D]. Jiangxi: East China University of 
Technology, 2022. 

22 Paul K. Romano, et al., OpenMC: a state-of-the-art Monte Carlo code for research and development, Ann. Nucl. Energy 82 (2015) 
90—97. 

23 WER, f) Jk BC 8] ZR ERES dE 
78,2022,45(01):73-81. 
Shen Zhirui, Sun Qizheng, He Donghao et al. Comparative verification of OpenMC program and NECP-X program based on high 
fidelity modeling of BEAVRS benchmark problem [J]. Nuclear Technology, 2002,45(01):73-81. 

24 吴 宏 春 . 中 子 输 运 理论 及 其 数值 方法 [M]. 西安 : 西安 交通 大 学 出 版 社 , 2007: 34. 
WU Hongchun. Neutron transport theory and numerical methods[M]. Xi'an: Xi'an Jiaotong University Press, 2007: 34. 

网 格 的 穿 透 概率 方法 研究 [J]. 核 动力 工程 , 2006, 27(5): 12-18, 41. DOI: 10.3969/j.issn.0258- 


IL 


P BEAVRS 基准 题 高 保 真 建 模 的 OpenMC 程序 和 NECP-X FEF X E SS uET]. Ez dx 


NS 


25 XE, RER. 基于 三 
0926.2006.05.003. 

LIU Pingping, WU Hongchun. Study of transmission probability method based on triangle meshes[J]. Nuclear Power Engineering, 
2006, 27(5): 12-18, 41. DOI: 10.3969/j.issn.0258-0926.2006.05.003.[ 知 网 中 文 ][ 知 网 英文 ] 

26 刘 国 明 ， 吴 宏 春 . 基于 三 棱柱 网 格 的 穿 透 概率 方法 研究 [中 . 核 动 力 工程 , 2008, 29(4): 24-30. 

LIU Guoming, WU Hongchun. Transmission probability method based on triangular-z mesh[J]. Nuclear Power Engineering, 2008, 
29(4): 24-30.[/3 77] 

27 BRE, 高 伟 ， 唐 富 初 , 等 . 三 维 中 子 积 分 输 运 方程 的 数值 解 : 中 子 碰撞 概率 的 计算 [中 .工程 数学 学 报 , 1984, 1(1): 119-128. 
MA Yichen, GAO Wei, TANG Fuchu, et al. Numerical solution of three-dimension integral transport equation: calculation of 
neutron collision probabilities[J]. Chinese Journal of Engineering Mathematics, 1984, 1(1): 119-128.[ 知 网 中 文 ][ 知 网 英文 ] 

28 Technical report : G. Marleau et al. AUSER GUIDE FOR DRAGON VERSIONS[B]. IGE-335, 2021. 


IIT 


29 Yamamoto A, Tabuchi M, Sugimura N, et al. Derivation of optimum polar angle quadrature set for the method of characteristics 
based on approximation error for the bickley function[J]. Journal of Nuclear Science and Technology, 2007, 44(2): 129—136. DOI: 
10.1080/18811248.2007.9711266.[LinkOut] 


XXXXXX-13 


